Quantum Tunneling in Half-Integer-Spin Kagome-Lattice Langasites 



A. Zorko,i'2 F. Bert,'"' P. Mendels,^'^ A. Potocnik,! A. Amato,^ C. Baines,^ 

K. Marty,s P. Bordet/ P. Lejay,"^ E. LhotelJ V. Simonet/ and R. Ballou^ 

^ Jozef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia 
^EN-FIST Centre of Excellence, Dunajska 156, SI-1000 Ljubljana, Slovenia 
^ Laboratoire de Physique des Solides, Umversite Paris-Sud 11, UMR CNRS 8502, 91405 Orsay, France 
^Institut Universitaire de France, 103 Boulevard Saint-Michel, F-75005 Paris, France 
"^Laboratory for Muon Spin Spectroscopy, Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 
^SPINTEC, (UMR8191 CEA/CNRS/UJF/G-INP), CEA Grenoble, INAC, 38054 Grenoble, France 
''institut Neel, GNRS and Universite Joseph Fourier, BP 166, 38042 Grenoble, France 

(Dated: November 1, 2012) 

Employing the muon spin relaxation technique we evidence temperature independent mag- 
netic fluctuations persisting down to the lowest temperatures in the samarium-based (J = 5/2) 
kagome-lattice Langasite. A detailed bulk-magnetization characterization and comparison to the 
neodymium-based (J = 9/2) compound allow us to assign the persistent spin dynamics to a quan- 
tum tunneling process. This is facilitated by pairwise anisotropic magnetic interactions, leading to 
a universal scaling of the muon relaxation. Our study reveals a remarkable analogy between weakly 
interacting half-integer-spin rare-earth magnets and molecular nanomagnets. 

PACS numbers: 76.75.+i, 75.45.+j, 75.10.Kt, 75.30.Gw 



Quantum tunneling (QT) mechanism allows transi- 
tions between different metastable ground states (GSs) 
through a potential barrier [l[ , which represents a salient 
alternative to classical thermally activated jumps over 
the barrier. It has been deeply investigated in molec- 
ular nanomagnets (MNMs) containing several magnetic 
ions entangled into a collective spin state 0-131, because 
of their high application potential in spintronics jlj. Re- 
cently, lanthanide compounds have been highlighted as 
mono-nuclear systems that may magnetically behave like 
MNMs d, 01 . Similarly to many MNMs they generally 
possess large magnetic moments and large inherent mag- 
netocrystalline anisotropy (compared to magnetic inter- 
actions) setting the potential barrier, which can even be 
tuned by chemical design to match challenging require- 
ments of quantum computation Q . 



by ^SR |12| and neutron spin-echo (NSE) experiments 
hj below 10 K, i.e., a temperature seemingly unrelated 
to at least an order of magnitude smaller exchange in- 
teraction and an order of magnitude larger crystal-field 
gap between the GS Kramers doublet and the first ex- 
cited state [H, 13|. The fluctuating GS was found to be 
surprisingly sensitive to the applied magnetic field, as a 
field-induced crossover from a dynamical to a more static 
state was observed |12| and even field-induced partial or 
short-range order was proposed 15, [l6| . Furthermore, 
zero-field low-temperature heat-capacity data suggested 
a Schottky anomaly that would require splitting of the 
GS doublet Despite all the efforts a comprehensive 
explanation of these unusual findings is still missing. 

In this Letter, we unveil that the persistent spin dy- 
namics are intrinsic to half-integer-spin Langasites as well 



In our study we have focused on the neodymium- 
and samarium-based members of the Langasite family 
R3BC3D2O14 [R - rare earth) 0, the first physical re- 
alizations of a geometrically frustrated kagome lattice 
(Fig. [1]) with dominant magnetocrystalline anisotropy. 
The Nd-based Langasite Nd3Ga5SiOi4iNGS) has been 
intensively studied in recent years [71, Il0l4la | due to its 
possible spin-liquid GS. Muon spin relaxation (/xSR) 
|12| alongside neutron scattering |ll| indeed unambigu- 
ously demonstrated the absence of electronic-spin freez- 
ing down to at least 60 mK. Since large antiferromag- 
netic exchange interactions between Nd'^"'" magnetic mo- 
ments (J = |) were initially predicted 0|, this was at- 
tributed to strong geometrical frustration. Subsequent 
refined investigations, however , suggested exchange cou- 
pling in a sub-Kelvin range [l3|, |l4| , although the exact 
value could not be determined. A plateau in spin relax- 
ation evidencing persistent spin dynamics was detected 




FIG. 1: Distorted kagome network of rare-earth R'^'^ magnetic 
moments in Langasites, R^BC^DiOiA [ab plane at c = 0). 
Contour shows the electrostatic potential of La3Ga5SiOi4 
with blue (magenta) regions corresponding to high (low) val- 
ues. Positions Pi and P2 are most favorable to muons. 
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as we provide further insight to the unusual magnetic- 
field dependence of the magnetic fluctuations. We fo- 
cus on the highly sensitive local-probe fj.SK technique 
Jl that has already proven invaluable in Langasites 
I2I . llSfl . Investigating the new Kramers- ion represen- 



tative Sm3Ga2.63Al2.37SiOi4 (SGS) with Sm^+ (J = |) 
magnetic moments along with combining new and previ- 
ously existing ^SR data on NGS enables us to clarify 
the origin of the persistent spin dynamics. These arise 
from the quantum tunneling process, leading to a uni- 
versal scaling of the field-dependent relaxation with the 
magnitude of the magnetic moments. 

In our investigation, the same NGS powder samples 
were used as previously [l^, while the SGS pow- 
der samples were synthesized for the first time [l^. Our 
syntheses led to highly crystalline samples with the dom- 
inant phase Sm3Ga2.63Al2.37SiOi4, and minor impurity 
phases Sm3Ga50i2 and Sm4.66 0(Si04)3 Bulk mag- 
netic properties between 1.6 and 830 K were checked 
with a commercial Quantum Design MPMS SQUID mag- 
netometer and two purpose-built axial extraction mag- 
netometers. Low-temperature measurements down to 
60 mK were performed with a home-built SQUID magne- 
tometer equipped with a miniature dilution refrigerator. 

Like in NGS, bulk susceptibility (x = AI/H in the 
linear regime; M is magnetization and H is the applied 
magnetic field) shows no indications of spin freezing in 
SGS down to 2 K [Fig.Eta)]. A maximum in (M/H)-^ 
observed around 400 K is due to the proximity of the first 
excited ^Ht and the GS ^Hs J-multiplet of Sm^+j^^. 
The latter is further split by a low-symmetry local crys- 
tal field. Therefore, at low temperatures only the GS 
Kramers doublet is populated, yielding an effective S — ^ 
spin state for both NGS and SGS. This is evidenced 
by the low-temperature field dependence of magnetiza- 
tion that fits nicely to the normalized Brillouin function 
^■Bi_{iJLgB/kBT) [¥ig.mh)l where fig and fij = gjfigJ 
are the GS and the free-ion magnetic moment, fj,B is 
the Bohr magneton, ks the Boltzmann constant and 
Bi{x) = tanh(a;). The fit yields magnetic moments 
^^NGS ^ i Q^^ j^jjd ^SGS ^ o.25(3)/iB. To obtain 
the latter we take into account powder averaging and as- 
sume analogy to NGS where the moment perpendicular 
to the c axis is ^^^^ = O.Sfi^'^^ 0. The suitability of 
the Brillouin function indicates weak magnetic interac- 
tions. In NGS, these can be evaluated by comparing the 
low-temperature susceptibility of the pure and a diluted 
NGS sample [inset in Fig. [^^a)], in which non-magnetic 
La'^"'" ions randomly substitute for 99% of the Nd'^^ ions 
The Curie- Weiss behavior {M/H)-^ = {T -9)/C 
yields the Weiss temperature 9 = —120(20) mK in the 
pure and ^ in the diluted sample. 

Invaluable local insight to magnetism in both materials 
is provided by /iSR. Measurements were performed on the 
General Purpose Surface-Muon (GPS) and the Low Tem- 
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FIG. 2: (a) Temperature dependence of {M/H)~^ (M is 
magnetization) in Sm3Ga2.63 Al2.37SiOi4 (SGS) measured in 
the field hqH — 1 T below 350 K and extracted from the 
M(H) curves above 350 K. Inset shows {M/H)~^ measured 
in 10 mT in Nd3Ga5SiOi4 (NGS) and magnetically diluted 
(Ndo.oiLao.gg)3Ga5SiOi4, with the linear extrapolation inter- 
secting the temperature axis, (b) Field dependence of magne- 
tization normaUzed by the saturated value for free ions mea- 
sured at 2 K on a powder SGS sample and at 1.6 K on a NGS 
single crystal (H||c) Solid lines are fits to the Brillouin 
function ^ ■ Bi{figB/kBT). 



perature Facility (LFT) instruments at the Swiss Muon 
Source (SfiS), Paul Scherrer Institute (PSI), Switzer- 
land. In fiSK experiments the spin of almost 100% spin- 
polarized muons is used to locally probe the magnetic 
field at a muon stopping site, causing precession of 
the muon spin. Static fields lead to coherent oscillations 
of muon polarization. Their distribution dampens the 
oscillations, leading to a single dip in the case of ran- 
dom frozen fields On the other hand, for fast fiuc- 
tuations of [v ^ where v is the fluctuation 

frequency and 7^ = 851.6 MHz/T is the muon gyromag- 
netic ratio) the muon polarization generally decays ex- 
ponentially, P = c^"^*, in zero (ZF) or longitudinal mag- 
netic field (LF) -Blf- The muon spin- lattice relaxation 
rate A = /p ^'^{B±{t)B±{Q))cos{jf^Bi,-pt)dt measures a 
spectral- width of spin fiuctuations through correlation of 
the local magnetic field B± perpendicular to the applied 
field. Description of these correlations with a single cor- 
relation time yields the Redfield's relation Q 



A 



(1) 



LF 



As previously shown for NGS [12|, also for SGS the 
experimental decay of the muon polarization can be fit- 
ted with the stretched-exponential relaxation function 
P{t) = e^^'^*^ with constant a = 0.6(1) in the range 
between 20 mK and 100 K [Fig. [21Ia)[, once the small 
nonrelaxing background part ^0.05 is subtracted. The 
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stretch exponent a < 1 is explained by multiple stopping 
sites in Langasites, as argued later on. The monotonic 
decay immediately discloses dynamical local fields even 
at the lowest experimentally accessible temperatures. T- 
dependence of A measured in ZF reveals strong varia- 
tions of spin fluctuations at high temperatures in both 
compounds. The observed activated behavior A oc 1/;^ oc 
e^/"^ [inset in Fig. [3] is expected for T < A, where A 
denotes the splitting of the lowest excited state from the 
GS doublet. It is due to crystal-field transitions driven 
by a magneto-elastic coupling Q|. This single-ion mecha- 
nism leads to vanishing relaxation in the T ^ limit. In 
contrast, we observe saturation of the muon relaxation 
[Fig- E^b)], which, quite interestingly, occurs in SGS be- 
low the same temperature Tc 10 K as previously found 
in NGS [l^ . The relaxation plateau evidenc es q uantum 
spin fluctuations beyond the classical picture [ISj setting 
up in both compounds and was observed before in several 
rare-earth-based frustrated magnets 21-2^. 

The entire T-dependence of A in both NGS and SGS 
can be fitted to the phenomenological model 2^ 2^ 



A(r) = 



A/T ' 



(2) 



valid for fast fluctuations when A oc l/v. In this 
model both the T-dependent phonon-induced broadening 
y4g-A/T ^^Y[ as the T-independent quantum contribu- 
tion 1/A„j are incorporated in the spectral width of spin 
fluctuations. Our flts [Fig. E^b)] yield the low-T relax- 
ation rates A^'^^ = 10.7 fJ-s^^, Xf^^ = 2.5 fis^^ and acti- 
vation gaps A^cis ^ K, A^gs ^ 2IO K for NGS and 
SGS, respectively. The dynamical energy barrier A'^^^ 
agrees with the NSE experiment and is consistent with 
the static crystal-field splitting (90 K) between the GS 
and the lowest-lying crystal- field doublet [l^. We find 
that in SGS this splitting is even larger. 

Next, we elucidate the magnetic GS of NGS and SGS 
and highlight the mechanism responsible for quantum re- 
laxation leading to the relaxation plateau. Our density- 
functional-theory (DFT) calculations predict interstitial 
sites Pi = (0.28,0.22,-0.08) and P2 = (0.27,^0.68,0) 
as the most probable muon stopping sites [19|, since 
electrostatic minima are found there (Fig. [T]). An es- 
timate of the fluctuation rate v = 2j'^B'j^/Xm [Eq- (HI] 

can be obtained by calculating the ensemble-averaged lo- 

1/2 

cal magnetic fleld B± = {B\) . Since Langasites are 
insulators with well localized 4/ electrons, the dipolar 
contribution dominates. For NGS our calculations, tak- 
ing into account noncorrelated moments /i^*^^ = 1.6/iB 
within a sphere with a radius large enough to ensure 
convergence, yield B^^ = 0.11 T and B^^ = 0.19 T 
for Pi and P2, respectively. With X^^^^ = 10.7 ns'^ 
these fields give on average j/^GS _ GHz. a sim- 
ilar value as in our previous estimate for oxygen sites 
[1^ . This is in excellent agreement with the NSE ex- 
periment, which determined the spin correlation time 




0.01 



2 4 



6 0.1 1 10 100 
7-(K) 



FIG. 3: (a) Relaxation of ZF muon polarization in 
SmsGaa.esAla.srSiOw (SGS) and NdgGasSiOw (NGS) [Hj 
at selected temperatures. Lines correspond to the stretched- 
exponential fit P(t) = e~'^*' . (b) T-dependence of the muon 
spin- lattice relaxation rate with fits (lines) to Eq. ((2)). Inset: 
the activated behavior A oc e^^"^ at high temperatures. 



^NGS _ ^ Q 3 jjg_ Besides putting confidence 

in our DFT calculations, this accordance of complemen- 
tary techniques confirms that muons directly probe the 
relaxation of electronic moments in Langasites and rules 
out any possible muon induced effect. The observation 
of the juSR relaxation plateau, therefore, proves intrinsic 
persistent spin dynamics in Langasites. The same anal- 
ysis in SGS, taking into account the above determined 
0.25(3)AiB and A^ps ^ 2.5 fis'^, yields the fluc- 



,,SGS 
^5 



tuation rate i/ 



SGS _ 
,NGS 



0.36(5) GHz, which is suppressed 
compared to 1/^"'= by fj.^'^^ / ^j,^^^ . 

Our previous /iSR investigation of NGS witnessed a 
pronounced fleld-induced maximum of A at 60 niK [l^ . 
contradicting Eq. ([T]) that would predict a monotonic de- 
crease of A if the fluctuating rate i' was constant. We find 
that the maximum at P^^s = 0.3 T is T-independent 
throughout the plateau region (Fig. [5]). Therefore, field 
polarization cannot play any role in its formation. Since 
7^i?^^^ <C v{Bi,F = 0) the initial increase of A with 
increasing Blf reveals decreasing fiuctuating rate that 
reduces by a factor of 2 at B^'^^ . This observation con- 
tradicts the traditional paramagnetic picture, in which a 
monotonic increase of ly due to a reduced energy barrier 
between the ground and excited states in the applied field 
would nionotonically decrease A. Similar maxima in 
A have been predicted and observed jUl in MNMs 
and should be explained within the energy-level scheme 
of the system. 

^NGS gg^g relevant energy scale in NGS, J'NGS _ 
^NGS^NGS^^^ « 300 mK, which nicely corroborates 
with inelastic neutron scattering [ll| and specific-heat 
measurements [l^j both suggesting ZF splitting of the 
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FIG. 4: Field dependence of the muon relaxation rate A in 
Sm3Ga2.63Al2.37SiOi4 (SGS) and NdgGasSiOw (NGS). Lines 
correspond to the model given by Eq. ([l}, with ZF parameters 
for NGS (dashed line) and SGS (solid line). Inset shows the 
universal scaling of A/^j with Bi^p/ng. 



GS Kramers doublet by ~250 mK. This suggests that 
j-NGS jg associated with magnetic anisotropy, because 
the isotropic exchange cannot split the Kramers dou- 
blet. Such magnetic anisotropy cannot be ascribed 
to the dipolar interaction that is too small, J^'^^ = 
^^'^^Bd/ks ~ 25 mK (average dipolar field is Bd = 
25 mT). This is further confirmed by the dipolar con- 
tribution to the Curie-Weiss temperature [29| 6d = 

1 ^ NGS / NGS^4 o 2 

l-(A'i /Mg ) 3z, 



= — 13 mK that 



is an order of magnitude below the experimental value 
6 = —120(20) mK. Thus, we attribute the energy scale 
j-NGS anisotropic part of the tensorial exchange 

interaction, which is for rare earths with large spin-orbit 
coupling often of the same magnitude as the isotropic 
part [30| . We note that 6 includes contributions of both 
parts and is, therefore, also in agreement with J'^*^^. 

In SGS we find a similar field dependence of the re- 
laxation rate, provided that both B^p and A are nor- 
malized by the moment value fXg (inset in Fig. U]). This 
evidences the characteristic scale Bc/^g ~ 0.1 T/fis, in- 
dependent of the rare-earth ion. It yields J'SGS^j-ngs ^ 

'-'^) , further supporting our claim that the en- 
ergy scale should be attributed to a pairwise inter- 
action. In this case, the expected dependence v/ng = 

^ leads to the scaling 

-PTTs — 1 — TTTs — -/ — \2 ■ The pairwise interaction 
i7 should bring about building up of spin correlation for 
T < J . This could indeed be responsible for the promi- 
nent decrease of }^'^^ below '-^0.5 K [Fig. El^b)], which 
is absent in SGS in the same temperature range because 
the relevant energy scale j/^SGS jg much smaller. 



/(Slf/Ms) according to Eq. 



Persistent fluctuations at temperatures T/A ^ 10~^ in 
NGS and SGS can be explained by QT through a double- 
well potential separating the two states of the crystal- 
field Kramers doublet. A matrix element that couples 
these two states and induces tunneling has to be asso- 
ciated with magnetic anisotropy, resulting in a splitting 
of the GS. In MNMs various sources of anisotropy have 
been invoked, such as transverse single-ion anisotropy, 
hyperfine interaction, dipolar fields and anisotropic ex- 
change interactions [2| . Since in NGS and SGS the muon- 
relaxation plateau persists to temperatures far above J ^ 
single-ion spins are tunneling due to quasi-static local 
magnetic fields that break the time-reversal symmetry. 
The single-ion anisotropy that leads to QT in integer-spin 
MNMs a is ineffective in half-inte ger-spin Langasites 
due to topological reasons jsil . Is^ . In half-integer-spin 
MNMs quasi-static nuclear hyperfine fields have proven 
responsible for muon relaxation (26l . |33| . However, these 
cannot explain QT in Langasites, because the hyperfine 
interaction in Nd'^+ and Sm'^"'" is of the same magnitude, 
A J ^ 10 mK [34|, which should lead to very similar fluc- 
tuation rates j/^^is ^^^^ j/SGS. g ^^g fQu^jj foj. a family 

of isotropic MNMs [s^- Moreover, in both materials the 
natural abundance of magnetic rare-earth nuclei is be- 
low 30%; thus the stretch exponent a would be different 
inside and outside of the QT regime. Finally, since the 
dipolar interaction is much smaller than , we assign 
the anisotropic exchange interaction between rare-earth 
moments as being responsible for QT in Langasites. It 
produces longitudinal as well as perpendicular effective 
magnetic fields, the latter inflicting QT (sst . 

We propose a picture of tunneling under "quenched" 
anisotropic-exchange fields, where stochastic fluctuations 
of the environment are slow compared to spin correlation 
times T^*^^ = 0.33 ns and r^^^ = 2.7 ns of reference 
spins. This condition is satisfied when the thermally- 
assisted spin relaxation becomes slower than the quan- 
tum relaxation [Eq. which occurs around 20 K in 
both systems and thus rationalizes the materialization of 
the quantum relaxation at temperatures two orders of 
magnitude above the dominant interaction J . 

In conclusion, we suggest that the persistent spin dy- 
namics in the Kramers representatives of the kagome- 
lattice Langasites are due to QT and are triggered by 
short-ranged magnetic anisotropy splitting the GS dou- 
blet. Similarly, a tunneling process between different 
equivalent spin configurations has been proposed to ex- 
plain the dynamics of spin-ice pyrochlores above freezing 
QT favors a superposition of degenerate states 
over selecting a particular state. Therefore, it may af- 
fect frustrated lattices profoundly, possibly leading to a 
quantum spin-liquid GS (38| . like the dynamical quan- 
tum spin-ice state 3^ expected to materialize in sev- 
eral pyrochlores 33^411 with substantial anisotropic ex- 
change interaction. QT thus appears as a momentous 
feature of the geometrically frustrated rare-earth-based 
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compounds and a close resemblance exists between these 
magnetically anisotropic systems and mesoscopic molec- 
ular nanomagnets. 
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Synthesis and Crystal Structure 

Powders of NGS and SGS compounds were prepared 
by solid state reactions of stoichiometric amounts of high 
purity oxides at 1420°C, resp. 1400°C in air. X-ray pow- 
der diffraction patterns were measured to confirm the 
crystal structures reported in Ref. In the SGS Lan- 
gasite with a heavier rare earth atom, part of Ga has 
to be replaced by the smaller Al in order to stabilize 
the structure. The Rietveld refined diffraction pattern 
of SGS is shown in Fig. [S] and the structural parame- 
ters are reported in Tab. ID Except from the Sm site, all 
other cation sites show substitutions. The Si cation oc- 
cupancy is fixed by charge balance, and the refinement of 
the Ga/Al mixed occupancies led to the overall compo- 
sition for SGS; Sm3Ga2.63Al2.37SiOi4. The NGS single 
crystal was grown by the floating zone method using an 
image furnace, under a 99% Ar + 1% O2 atmosphere, at 
a growth rate of 10 mm/h. 



Determination of Muon Stopping Sites 

Positively charged muons favor electronegative envi- 
ronment. In oxydes, it is generally believed that a muon 
stops close to oxygen, at a distance of about 1 A In 
order to determine possible muon stopping sites in Lang- 
asites we employed density-functional-theory (DFT) cal- 
culations for the isostructural La3Ga5SiOi4 compound. 
The goal was to determine a self-consistent electron- 
density distribution, which yielded a spatial profile of the 
electrostatic potential. Similar calculations have proven 
successful in the past 0|. 

Our calculations were performed using the pwscf pro- 
gram, a part of the Quantum Espresso software pack- 
age [6|. Ultrasoft pseudo potentials appropriate for 
the Perdew-Burke-Erzerhof exchange-correlation (LDA) 
were used. The pseudo Bloch functions were expanded 
over plane waves with the energy cutoff of 35 Ry on a 
6x6x6 Monkhorst-Pack fc-space mesh. The convergence 



3600 r- 
3100 - 
2600 - 




260 



FIG. 5: SGS x-ray diffraction pattern at room temperature 
(red dots) with its Rietveld refinement (black line) showing 
the presence of 88.7(6) mass % of the Sm3Ga2.63 Al2.37SiOi4 
Langasite, 5.8(1) mass % of the Sm3Ga50i2 garnet and 5.5(1) 
mass % of Sm4. 660(8104)3 01, as well as some contamination 
of the Cu sample holder. The difference between the mea- 
surement and the calculation is shown as a blue line. 

of the results and the total energy was tested and yielded 
negligible variation at this values. The Fermi surface was 
smeared with a "temperature" parameter of 2 meV. The 
computed electrostatic potential consisted of bare ionic 
and Hartree electrostatic potentials. 

A global electrostatic potential minimum was found at 
the interstitial site Pi = (0.28, 0.22, -0.08) with multi- 
plicity of six (6g), while local electrostatic minima were 
found at P2 = (0.27,0.68,0) and P2' = (0.59,0.32,0), 
both with multiplicity of three (see Fig. 1 of the main 
part). The value of the potential at P2 is only slightly 
enhanced with respect to the global minimum. Differ- 
ent value of the potential at crystallographically equiv- 
alent P2 and P2' is due to their proximity to the Ga(3) 
crystallographic site, which is randomly occupied with 
Ga?+ and Si4+ 0. In the case of Si^"*" the local mini- 
mum becomes significantly less intense. As expected, all 
electrostatic minima are located in a vicinity of the most 
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TABLE I: Structural parameters of SGS determined by x-ray powder diffraction with a D8 Brucker diffractometer. The data 
refined with the Rietveld method using the Fullprof software Q led to Bragg R-factor= fO.9, Rf-factor= 8.35, x^ = l-4 and the 
cell parameters a = 7.94798 (3) A, c = 4.96578 (3) A at room temperature. The Ga3/Si3 site is half occupied by Si and a 
mixture of Ga/Al. The columns are the atom type, the Wyckoff site, the atomic positions, the isotropic displacement parameter 
and the occupancy. 



Atom 


Site 


X 


Y 


Z 




Occ. 


Sm 


3e 


0.4157(2) 








0.5 


1 


Gal 


la 











1.074(7) 


0.432 


All 


la 











1.074(7) 


0.564 


Ga2 


3f 


0.7706(9) 





1/2 


1.074(7) 


0.536 


A12 


3f 


0.7706(9) 





1/2 


1.074(7) 


0.464 


GaS 


2d 


1/3 


2/3 


0.4603(31) 


1.074(7) 


0.294 


A13 


2d 


1/3 


2/3 


0.4603(31) 


1.074(7) 


0.204 


Si3 


2d 


1/3 


2/3 


0.4603(31) 


1.074(7) 


0.501 


01 


2d 


1/3 


2/3 


0.8025(54) 


1.0 


1 


02 


6g 


0.4638(40) 


0.3169(23) 


0.7032(28) 


1.0 


1 


03 


6g 


0.2301(22) 


0.0783(18) 


0.2514(22) 


1.0 


1 



electronegative ions. Pi at the distance of 1.27 A and 
P2 (-Fb') 0.93 A away, in agreement with the generally ac- 
cepted distance Assuming that the rather localized 
4/ electrons do not appreciably affect the electrostatic 
picture except in the close vicinity of rare earths and that 
the muons do not perturb the crystal structure, we assign 
Pi and P2 as the most probable muon stopping sites also 
in Nd3Ga5SiOi4 and Sm3Ga2.63Al2. 3781014. The same 
stretched-exponential functional form of the muon depo- 
larization curves (a = 0.6; see main part) in both ma- 
terials suggests that small structural differences between 
them have no major impact on the muon stopping sites. 
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